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Abstract 



Lattice QCD with Wilson fermions generically shows the phenomenon of a first order 
phase transition. We study the phase structure of lattice QCD using Wilson twisted 
mass fermions and the Wilson plaquette gauge action in a range of (3 values where such 
a first order phase transition is observed. In particular, we investigate the dependence 
of the first order phase transition on the value of the lattice spacing. Using only data 
in one phase and neglecting possible problems arising from the phase transition we are 
able to perform a first scaling test for physical quantities using this action. 



1 Introduction 



Understanding the phase structure of lattice QCD is an important pre-requisite before start- 
ing large scale simulations. Indeed, our collaboration found that when working at lattice 
spacings of about 0.15 fm there can be strong first order phase transitions at small quark 
masses, at least when a combination of Wilson plaquette action and Wilson fermions is used 
PIE]- The phenomenon appears also when a small twisted mass term is switched on. This 
has serious consequences, since in such a scenario the pion mass mps cannot be made arbi- 
trarily small but assumes a minimal value, m™, which may be about 500 MeV and hence 
it becomes impossible to work close to the physical value of the pion mass. 

The presence of the first order phase transition for pure Wilson fermions is in accordance 
with predictions from chiral perturbation theory [Sj, which have been extended later to 
the case of adding a twisted mass |U El El El E]- Let us, for completeness, also mention 
that for values of the lattice spacing much coarser than a = 0.15 fm the first order phase 
transition turns into a second order one from the normal QCD phase to the so-called Aoki 
phase |H1 EH1 EI]. The generic phase structure of lattice QCD according to our present 
understanding is discussed and illustrated in refs. [Tl 121 IT2]. 

In refs. [U E] we have studied only one value of the inverse gauge coupling (3 = in 
order to demonstrate the existence of the first order phase transition, leaving the question 
of the (3 dependence open. Since lattice chiral perturbation theory predicts a weakening of 
the first order phase transition towards the continuum limit, it is interesting to check this 
prediction and, in particular, to investigate quantitatively how fast the transition weakens 
when the continuum limit is approached. The answer to the latter question will naturally 
depend on the choice of the actions that are used for the gauge and the fermion fields. 

In this paper we will present results using Wilson twisted mass fermions and the Wilson 
plaquette gauge action for three values of (3. At each of these (3 values we have performed 
simulations at a number of quark masses on both sides of the first order phase transition. 
This allows to study the (3 dependence of the phase transition itself and, in addition, the 
lattice spacing dependence of physical observables computed separately in the two phases. 
We have performed such a scaling test for the pion mass, the pion decay constant and the 
ratio of the pion to the vector meson mass. For a scaling test of Wilson twisted mass fermions 
and other recent results in the quenched approximation see refs. [131 EU ES] . 
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2 Wilson twisted mass fermions 



In this paper we will work with Wilson twisted mass fermions [M] that can be arranged to 
be 0(a) improved without employing specific improvement terms [Tj\. The Wilson tmQCD 
action in the twisted basis can be written as 

S[U,x,x] = a 4 ^2x(x)(D w + m + i/j,>y 5 T 3 )x(x) , (1) 

X 

where the Wilson-Dirac operator Dw is given by 

3 1 

= 2 + V ^ " aV M V ^ (2) 

At=0 

and V M and V* denote the usual forward and backward derivatives and the Wilson parameter 
r was set to 1. 

The situation of full twist and hence automatic 0(a) improvement arises when mo in 
eq. (UJ) is tuned towards a critical bare quark mass m crit . We use for our simulations the 
hopping representation of the Wilson-Dirac operator with k = (2am + 8) _1 . 

We extract the pseudo scalar mass mps and the vector meson mass my from the usual 
correlation functions: 

C PP (x ) = a 3 J2( P+ ( x ) P ~(°)) > 

3 " (3) 
CM*o) = y]C]r<y+(*)V7(0)>, 

fc=l X 

where we consider the local bilinears P ± = xis^X an d = Xln^X- Here we used 
= (ti ±zt2) with 7~i ; 2 the first two Pauli matrices. Similarly one can define the correlation 

function Cap with the local bilinear = xipTlfr^X- 

The bare pseudo scalar decay constant in the twisted basis can be obtained from (cf. 

da CHI) 

f? = m£r AP (0\P + (0)\*), (4) 

where the ratio 

_ (0\At(0)\n) 
rAP ~ <0|P+(0)|tt) (5j 



2 



p 


L 3 x T 


a/i 


a [fm] 


5.1 


12 3 x 24 


0.013 


0.200(2) 


5.2 


12 3 x 24 


0.010 


0.160(4) 


5.3 


16 3 x 32 


0.008 


0.138(8) 



Table 1: Simulation points for Wilson plaqucttc gauge action. For the three values of /3 we give the lattice 
extent, the value for an and the value of the lattice spacing in fm, determined using ro = 0.5 fm at the 
reference point (see text), where (romps) 2 = 1-5. 



can be extracted from the asymptotic behavior of 

^4^4 = r AP tanh[m PS (T/2 - x )} . (6) 
The bare PCAC quark mass m^ CAC in the twisted basis can then be computed from the ratio 



fPS 

m * "2(0|P+(0)|7r) mps - (7) 



PCAC _ Jx 



The sign of m^ CAC and f^ s is determined by the sign of r^p and therefore, the corresponding 
values can be negative. One has to keep in mind that m PCAC and f^ s , since measured in 
the twisted basis, do not correspond to the physical quark mass and the physical pseudo 
scalar decay constant, respectively. While the quark mass is given by a combination of the 
(renormalized) values of m^ CAC and fi, the pseudo scalar decay constant can be computed 
by the help of f^ s and the twist angle, as long as f^ s ^ and the value of the twist angle is 
different from tt/2. 

Note that the purpose of the present paper is not to work at full twist nor to extract 
physical quantities, but rather to study the lattice spacing dependence of the first order 
phase transition. For the same reason, we also do not address the question of the choice of 
the critical quark mass in order to stay at full twist here, see refs. [TBI IT4] for recent quenched 
simulations addressing this point. 
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3 The phase transition as a function of the lattice spac- 
ing 



In order to study the lattice spacing dependence of the phase transition we have chosen three 
values of/3: /3 = 5.1, /5 = 5.2 and (3 = 5.3. We scaled the volumes and the values of /i such 
that the physical volume is larger than 2 fm, roughly constant and that r /i ~ 0.03, where 
tq is the Sommer scale [20] fixed to be tq = 0.5 fm. Note that the value of r^ja depends on 
the value of the quark mass and therefore we had to choose a reference value for r^/a as will 
be explained below. The parameters are summarized in table ^ 

In practice it turned out that a very direct way of detecting the presence of a first order 
phase transition in lattice QCD is to monitor the behavior of the plaquette expectation value 
(P), e.g. as a function of n for fixed twisted mass parameter /i. In such a situation, starting 
at identical parameter values from "hot" (random) or "cold" (ordered) configurations, (P) 
can assume different, co-existing values. In fig. Owe show (P) as a function of 1/(2k) for 
the three values of j3. The picture is typical for the behavior of a first order phase transition 
with meta-stable branches, one with a low value of (P) and one with a high value of (P). We 
will denote in the following these branches as high ("H") and low ("L") plaquette phases, 
respectively. 

The /3-dependence shows that the gap in the plaquette expectation value AP decreases 
substantially when moving from — 5.1 (a ~ 0.20 fm) to — 5.3 (a ~ 0.12 fm), which 
is presumably due to the mixing with the chiral condensate as discussed in £Q. One possi- 
ble definition for the quantity AP is the difference between low and high phase plaquette 
expectation value at the smallest value of k where a meta-stability occurs. 

Let us remark that the first order phase transition exists also in the continuum limit at 
zero quark mass where the scalar condensate has a jump as a consequence of spontaneous 
chiral symmetry breaking. This means, of course, that in the continuum limit the phase 
transition occurs only for fi = 0. 

We give our simulation parameters, the statistics of the Monte Carlo runs and the results 
for amps, a /^ S ' am x CAC an< ^ r o/ a i n tables EJ El and 

The meta-stability phenomenon observed in (P) can also be seen in fermionic quantities. 
As an example, we show in fig. |2] the values of the PCAC quark mass as obtained in the 
branches with high and low plaquette expectation values of fig. ^ for the three values of 
(3. Again we observe that with increasing /3 the gap between positive (low plaquette phase) 
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0.56 r ■ 




0.46 1 

2.80 2.84 2.88 2.92 2.96 3.00 3.04 



Figure 1: The plaquette expectation value (P) as a function of 1/(2k) at the three values of (3 we have 
simulated. We also indicate the values of a/i which are scaled with /3 such that tqh is roughly constant. The 
lines just connect the data points and only serve to guide the eye. For this study we have more simulations 
points than in the next figures and than in tables 01 Eland El 



and negative (high plaquette phase) quark masses shrinks. Also, the met a- stability region in 
l/(2/t) gets much narrower with increasing (5. 

The effects of the first order phase transition can also be seen in the pion mass and the 
value of the force parameter r$. We plot in fig. 01 an example of the pion mass as a function 
of the PCAC quark mass at (3 = 5.3. The most intriguing observation here is that due to 
the presence of the first order phase transition, the pion mass, say for positive quark masses, 
does not go to zero but rather reaches a minimal value, and jumps then to the phase with 
negative quark mass. This is, of course, just another manifestation of the jump in the PCAC 
quark mass in fig. |2I 

In fig. 01 we also show the values of r^/a in the low and high plaquette phases at (3 = 5.3. 
Note that the values of r Q /a are quite different when determined in the low and the high 
plaquette phases, which is a generic feature also for other values of j3 and even for different 
gauge actions, see ref. [T2j . 

An interesting question is, at which value of the lattice spacing a the minimal pion 
mass mpg n assumes a value of, say, 300 MeV where contact to chiral perturbation could be 
established. 

The pion mass assumes two different values for a fixed quark mass, once this quark mass 
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Figure 2: In the graph on the left the PCAC quark mass is plotted as a function of 1/(2k) at the three 
values of (3 we have simulated. Positive values correspond to the low plaquette phase while negative values 
correspond to the high plaquette phase. The statistical errors are on this scale for most of the points smaller 
than the symbols. In the right plot we give a closeup of the (3 = 5.3 results. 






m™ n [MeV] 


AP 


5.1 


>600 


0.0399(1) 


5.2 


>630 


0.0261(1) 


5.3 


>470 


0.0077(4) 



Table 2: Minimal pion mass mpg n in physical units in the low plaquette phase and AP for the three (3 
values. To set the scale we used tq — 0.5 fm and the value of ro/a measured for the corresponding simulation 
point. 



lies inside the meta-stability region. These two values for the pion masses correspond to the 
two phases that for a certain interval of quark masses co-exist. The precise determination 
of the meta-stability region is, of course, very difficult. We can, however, give an interval in 
k, [k!,k 2 ], that can be read from tables |4] and El for the three different values, where 
meta-stabilities occur in our simulation. In the following, we will mainly concentrate on the 
low plaquette phase since this is the natural choice for studying lattice QCD. Being interested 
only in the low plaquette phase we determine then a lower bound for the minimal pion mass 
as computed at the lower end of this interval, i.e. k\, in the low plaquette phase. We give in 
table 121 the values of the minimal pion masses in the low plaquette phase in physical units. 
In addition, we provide the value for the gap in the plaquette expectation value AP. 
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Figure 3: The squared pion mass as a function of the PC AC quark mass at j3 = 5.3. 



In principle, it would be interesting to extrapolate the minimal pion mass and the gap 
in (P) as a function of the lattice spacing. However, our present data do not allow for a 
reliable and safe extrapolation. First of all, the determination of the minimal pion mass has 
a large ambiguity in itself since we do not know exactly for which value of the quark mass the 
meta-stability will disappear. A substantially larger statistics would be necessary to answer 
this question and to check whether tunneling from one phase to the other occurs. Second, 
the only three values of (3 we have used give a too short lever arm to perform a trustworthy 
extrapolation. And, third, the values of r^/a are very different in the two phases, as can be 
seen in fig. 01 which makes it particularly difficult to follow the gap in (P) as a function of 
a/r . 

Nevertheless, an estimate on a more qualitative level yields a value of the lattice spacing 
of a ~ 0.07 fm — 0.1 fm where simulations with pion masses of about 300 MeV can be 
performed without being affected by the first order phase transition. 
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Figure 4: ro/a as a function of l/(2re) at ft = 5.3. 



2.993 



4 Lattice spacing dependence of 
physical observables 

Although the present simulations are not at full twist, the fact that we have results at three 
values of (3 with roughly constant r /i allows us to check for the size of lattice artifacts. 
In order to perform such an investigation it is advantageous to express physical quantities 
in dimensionless variables. To this end, let us first define a reference pion mass through 
(romps) 2 = 1.5. We have chosen this particular value in order to be able to interpolate for 
the values of (3 = 5.1 and (3 = 5.3, and to perform only a short extrapolation for (3 = 5.2 to 
this point. 

At the aforementioned reference pion mass, a corresponding reference value of r^/a and 
a reference quark mass can be determined, the latter leading to a variable a, 



m 



a 



PCAC 
X 



PCAC 



(8) 



I ref 



8 



p 


amZ CAC , 

X rcf 


af PS 

aJ X rcf 


( r o/a) rcf 


5.1 


0.035(2) 


0.195(06) 


2.497(29) 


5.2 


0.025(4) 


0.139(15) 


3.124(85) 


5.3 


0.022(1) 


0.122(07) 


3.628(60) 



Table 3: Reference values for am^ CAC , a/^ S and ro/a. The reference point is chosen such that (romps) 5 
1.5. The errors include the interpolation errors. 



Similarly, we can define ratios for a quantity O, 



R o = -^~ (9) 

I rcf 



where 0| re f is the quantity as determined at the reference pion mass. The values for several 
quantities at the reference point can be found in table 01 

In order to determine the reference values for m PCAC , f^ s and r , in a first step we 
interpolated m PCAC linearly as a function of (r mps) 2 to the point where (romps) 2 = 1.5 and 

extracted the reference value for mZ CAC ■ Then we determined the reference values for f PS 

x J x 

and ro by quadratically interpolating the data as a function of m PCAC to the reference value 
of m PCAC . We repeated the latter step with a linear interpolation finding agreement within 
the errors. The fits to the data have been performed with the ROOT and MINUIT packages 
from CERN (cf. [2B 122]), taking the errors on both axis into account. We remark that for 
the quantity r mpg we have neglected the correlation of the data between r /a and amps- 

For a given observable O, Ro is a universal function of a for fixed value of \i in physical 
units that allows for a direct comparison of results obtained at different values of (3 and, in 
principle, even for different actions. Deviations of results at different (3 values provide then a 
direct measure of scaling violations. In fig. |5]we show -R m 2 g as a function of a. Note that for 
the scaling analysis we take the data in the low plaquette phase only since this corresponds 
to the standard lattice QCD situation. We also remark that some of the points taken in this 
analysis might be meta-stable. Nevertheless, we assume here that these data can serve for 
checking scaling violations. Besides the data from the present work, we added also results 
from simulations at /3 — 5.6 [23], which were obtained, however, at vanishing twisted mass 
parameter \i = 0. 

A rather amazing consequence of fig. |3]is that, despite the fact that we are using coarse 
lattices, we cannot detect any scaling violation, at least within the (large) statistical errors 
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Figure 5: The dependence of the pion mass on the quark mass, i.e. R m ^ s as a function of a. Besides the 
data of this work, we added in the plot also results from Wilson fermion simulations at j3 = 5.6 which 
were obtained at jj, = 0. 



of our data. Even more, the results of our present simulations at small values of (3 agree 
with results from simulations with pure Wilson fermions at — 5.6 setting \x = 0. The 
same observation is made for -R/ps, see fig. E] and the ratio mps/my, see fig. These results 
indicate that the lattice artifacts and the effect of a non-vanishing twisted mass parameter 
fi are surprisingly small. We remark here that in the case of the ratios like Rm^ s an d Rf% s 
one could have cancellation of mass independent cutoff effects. One has also to have in mind 
that, due to the presence of the first order phase transition, the simulated pion masses are 
still larger than 500 MeV. Whether our findings also hold when one is approaching the chiral 
limit is certainly an interesting but open question. However, in a set-up with Wilson twisted 
mass fermions and Wilson plaquette gauge action this question cannot be answered at these 
values of the lattice spacing. 
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Figure 6: The ratio Rfps for the pseudo scalar decay constant as a function of a. We also added results 
from Wilson fermion simulations for Rf PS at f3 = 5.6 |23j obtained with /i = 0. 

5 Conclusions 



In this paper we have investigated dynamical Wilson twisted mass fermions employing 
the Wilson plaquette gauge action. We have performed simulations at three values of 
(3 = 5.1,5.2,5.3, corresponding to values of the lattice spacing of a w 0.20,0.16,0.14 fm, 
respectively. The non-zero values of the twisted mass parameter \x were chosen such that 
royU ~ 0.03 for all of the three (3 values. At these rather coarse lattice spacings we find clear 
signals of first order phase transitions that manifest themselves in a meta-stable behavior of 
the plaquette expectation value and fermionic quantities, such as the PCAC quark mass and 
the pion mass. 

We clearly observe that the gaps in quantities sensitive to the phase transition, such as the 
plaquette expectation value and the PCAC quark mass decrease substantially when f3 is in- 
creased. Unfortunately, with our present set of simulations, we are not able to quantitatively 
locate the value of the lattice spacing, where the effects of the first order phase transition 
becomes negligible and where a minimal pion mass of, say, 300 MeV can be reached. As 
an estimate of such a value of the lattice spacing we give a range of a ~ 0.07 fm — 0.1 fm. 
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Figure 7: The ratio mps/my as a function of a. Again, we also added results from Wilson fermion 
simulations at j3 = 5.6 |28| . 



Of course, this would mean that a continuum extrapolation of physical results obtained on 
lattices with linear extent of at least L = 2 fm would be very demanding, since the starting 
point for such simulations would already require large lattices. It is therefore very impor- 
tant to find alternative actions such that the value of the lattice spacing can be lowered 
without running into problems with the first order phase transition. One candidate for such 
an action, the DBW2 gauge action, is discussed in ref. [T2| where it has indeed been found 
that modifying the gauge action alone can substantially reduce the strength of the first order 
phase transition. We are presently investigating another possibility, the tree-level Symanzik 
improved gauge action j2U • 

Despite the problems arising from the presence of a first order phase transition, we per- 
formed a scaling analysis for the pion mass, the pion decay constant and the ratio mpg/my 
for the data obtained at the three values of (3 where we performed simulations. To this end, 
we only analyzed data from the low plaquette phase, since this is the natural choice for QCD 
simulations. 

By defining a reference pion mass at (romps) 2 = 1.5, we computed the ratio of mps and 
to the corresponding reference values as a function of the PCAC quark mass, again 
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measured with respect to the corresponding reference quark mass. We find that for these 
ratios the scaling violations are remarkably small and cannot be detected with the present 
precision of our data. Even more, when adding data from simulations of Wilson fermions 
with \i — at /3 = 5.6, then these data fall on the same scaling curve as our results on much 
coarser lattices and with twisted mass parameter switched on. This indicates that not only 
the lattice artifacts but also the effect of switching on a twisted mass of the order of r /x ~ 0.03 
are small, at least for the rather large pion masses simulated here. This finding is surprising 
since it suggests that continuum values of physical quantities can be already estimated from 
simulations at not too small lattice spacings. Of course, our scaling results suffer from the 
fact that they are obtained using data that might be meta-stable as a consequence of the 
presence of the first order phase transition. Hence, a scaling test with an action that does 
not lead to significant effects of the first order phase transition is mandatory to check the 
results presented in this paper. 
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1 v meas 




af PS 




Ta / O 


0.1758 


L 


160 


0.7015(031) 


+0.2856(60) 


+0.0799(12) 


2.178(8)(4)(20) 


U.l i Oo 


T 


1 fin 
loll 


U.bloo(U4U J 


+U.zooo(ob J 


+U.Uoy ^l^J 


Z.Zdo{o ){0){o) 


0.1765 


L 


160 


0.5353(068) 


+0.2201(76) 


+0.0446(16) 


2.370(12)(4)(26) 


0.1768 


L 


160 


0.4468(051) 


+0.1683(82) 


+0.0268(13) 


2.625(19)(22)(1) 


0.1758 


H 


160 


0.5323(126) 


-0.2065(119) 


-0.0496(25) 


3.926(26)(12)(10) 


0.1763 


H 


160 


0.6771(116) 


-0.2351(227) 


-0.0777(50) 


4.087(56)(4)(0) 


0.1765 


H 


160 


0.7231(111) 


-0.2595(232) 


-0.0864(26) 


4.053(18)(17)(3) 


0.1768 


H 


160 


0.7377(119) 


-0.2302(136) 


-0.0926(38) 


4.139(35) (16) (2) 


0.1770 


H 


160 


0.7530(189) 


-0.2212(189) 


-0.0977(59) 


4.045(28)(10)(4) 



Table 4: Parameters and physical observables for the simulations with /3 = 5.1. The lattice size in these 
runs was set to 12 3 x 24 and the twisted mass parameter to a/i = 0.013. We give the values for k and the 
number of measurements N mcas performed. We indicate with "L" or "H" whether the plaquette expectation 
assumes a low or a high value. Moreover, we give the values for mps, /^ S ' m x CAG an d r o in lattice units. For 
ro we give in addition to the statistical error two systematic errors, the first of them coming from possible 
excited state contaminations and the second from the necessary interpolation of the force in r. 



H 






amps 


< S 


am^ CAC 


r /a 


0.17125 


L 


320 


0.6057(025) 


+0.2289(35) 


+0.0650(08) 


2.618(20)(5)(49) 


0.17150 


L 


459 


0.5066(050) 


+0.1968(38) 


+0.0452(08) 


2.800(17)(9)(4) 


0.17175 


L 


320 


0.4189(071) 


+0.1540(84) 


+0.0292(17) 


3.038(28)(14)(4) 


0.17125 


H 


320 


0.4173(111) 


-0.1571(166) 


-0.0352(43) 


4.796(63)(65)(15) 


0.17150 


H 


318 


0.4220(126) 


-0.1566(219) 


-0.0349(50) 


4.282(61)(16)(0) 


0.17175 


H 


320 


0.4985(088) 


-0.1770(119) 


-0.0494(28) 


4.418(23)(23)(0) 


0.17250 


H 


320 


0.6462(131) 


-0.1974(087) 


-0.0874(24) 


4.767(51)(7)(3) 



Table 5: Parameter and physical observables for the simulations with (3 = 5.2. The lattice size in these runs 
was set to 12 3 x 24 and the twisted mass parameter to a^t = 0.01. See table^for further explanations. 
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H 




AT 

* » meas 


amps 


„ f PS 

a fx 


amf, CAC 

A, 


nr. 1 „ 

r /a 


0.16715 


L 


100 


0.3525(061) 


+0.1349(111) 


+0.0255(18) 


3.668(34)(9)(8) 


0.16720 


L 


101 


0.3460(075) 


+0.1259(063) 


+0.0233(15) 


3.551(47)(2)(1) 


0.16725 


L 


180 


0.3213(088) 


+0.1211(078) 


+0.0200(18) 


3.716(49)(24)(0) 


0.16730 


L 


243 


0.3208(037) 


+0.1160(063) 


+0.0191(15) 


3.730(35)(6)(9) 


0.16735 


L 


160 


0.2757(040) 


+0.0887(084) 


+0.0157(12) 


3.841(47)(28)(1) 


0.16730 


H 


388 


0.2656(054) 


-0.1037(093) 


-0.0140(12) 


4.84(10)(2)(1) 


0.16735 


H 


100 


0.3114(106) 


-0.1291(128) 


-0.0204(17) 


4.808(95)(33)(3) 


0.16740 


H 


100 


0.3027(094) 


-0.1237(137) 


-0.0210(19) 


4.703(90)(0)(1) 



Table 6: Parameter and physical observables for the simulations with j3 = 5.3. The lattice size in these runs 
was set to 16 3 x 32 and the twisted mass parameter to a/i = 0.008. For further explanations see table 0] 



References 

[1] F. Farchioni et al, Eur. Phys. J. C39, 421 (2005), [hep-lat/0406039]. 

[2] F. Farchioni et al, Nucl. Phys. Proc. Suppl. 140, 240 (2005), [hep-lat/0409098] . 

[3] S. R. Sharpe and R. Singleton Jr., Phys. Rev. D58, 074501 (1998), [hep-lat/9804028]. 

[4] G. Minister, JHEP 09, 035 (2004), [hep-lat/0407006]. 

[5] L. Scorzato, Eur. Phys. J. C37, 445 (2004), [hep-lat/0407023]. 

[6] S. R. Sharpe and J. M. S. Wu, Phys. Rev. D71, 074501 (2005), [hep-lat/0411021]. 

[7] S. R. Sharpe and J. M. S. Wu, Phys. Rev. D70, 094029 (2004), [hep-lat/0407025]. 

[8] S. Aoki and O. Bar, Phys. Rev. D70, 116011 (2004), [hep-lat/0409006] . 

[9] S. Aoki, Phys. Rev. D30, 2653 (1984). 

[10] E.-M. Ilgenfritz, W. Kerler, M. Miiller-Preussker, A. Sternbeck and H. Stiiben, Phys. 
Rev. D69, 074511 (2004), [hep-lat/0309057]. 

[11] A. Sternbeck, E.-M. Ilgenfritz, W. Kerler, M. Miiller-Preussker and H. Stiiben, Nucl. 
Phys. Proc. Suppl. 129, 898 (2004), [hep-lat/0309059]. 



15 



[12] F. Farchioni et al, accepted for publication in Eur. Phys. J. C (2005), [hep-lat/0410031], 
DOI 10.1140/epjc/s2005-02262-5. 

[13] X]F , K. Jansen, A. Shindler, C. Urbach and I. Wetzorke, Phys. Lett. B586, 432 (2004), 
[hep-lat/0312013]. 

[14] A. M. Abdel-Rehim, R. Lewis and R. M. Woloshyn, Phys. Rev. D71, 094505 (2005), 
[hep-lat/0503007]. 

[15] %l F j K. Jansen, M. Papinutto, A. Shindler, C. Urbach and I. Wet- 
zorke, accepted for publication in Phys. Lett. B (2005), [hep-lat/0503031], 
DOI 10.1016/j.physletb.2005.05.051. 

[16] ALPHA, R. Frezzotti, P. A. Grassi, S. Sint and P. Weisz, JHEP 08, 058 (2001), [hep- 
lat/0101001]. 

[17] R. Frezzotti and G. C. Rossi, JHEP 08, 007 (2004), [hep-lat/0306014]. 

[18] UKQCD, R. M. Baxter et al, Phys. Rev. D49, 1594 (1994), [hep-lat/9308020]. 

[19] F. Farchioni, C. Gebert, I. Montvay and L. Scorzato, Eur. Phys. J. C26, 237 (2002), 
[hep-lat/0206008]. 

[20] R. Sommer, Nucl. Phys. B411, 839 (1994), [hep-lat/93 10022]. 

[21] The ROOT system home page, http:/ /root.cern.ch/. 

[22] MINUIT home page, 

http: / / seal. web. cern.ch/seal/snapshot / work-packages / mathlibs / minuit /home. html. 

[23] C. Urbach, K. Jansen, A. Shindler and U. Wenger, hep-lat/0506011. 

[24] P. Weisz, Nucl. Phys. B212, 1 (1983). 

[25] P. Weisz and R. Wohlert, Nucl. Phys. B236, 397 (1984). 



16 



